On the two-photon decay width of the sigma meson 
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We shortly report on the two-photon decay width of the light <r-meson interpreted as a quarkonium 
state. Results are given in dependence on the cr-mass and the constituent mass of the light quark. 
The triangle quark-loop diagram, responsible for the two-photon transition, is carefully evaluated: a 
term in the transition amplitude, often omitted in literature, results in destructive interference with 
the leading term. As a result we show that the two-photon decay width of the a in the quarkonium 
$H ' picture is less than 1 keV for the physical range of parameters. 
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I. INTRODUCTION 

The two-photon decay of scalar mesons represents a valuable mechanism to possibly pin down their internal structure 
O ■ (see [J B, 0, B B B B B EH EH EH El) EH El] an d references therein). In particular, the transition of the scalar- 
isoscalar resonance a = / (600) into 77 has received much attention in the literature. It is commonly believed that a 
i-C ■ decay width of about 3-5 keV would favor a quarkonium interpretation of the /o(600). In this short work we aim to 
show that this is not the case: the decay width of a scalar-isoscalar quark-antiquark state, with flavor wave-function 

nn = J ^(uu + dd) and a mass between 0.4 and 0.8 GeV as favored by recent studies (a mass of about 0.44 GeV is the 

■ outcome of [IB]), turns out to be smaller than 1 keV for the physical range of parameters. When evaluating the related 
quark triangle-loop diagram of Fig. 1 care has to be taken concerning gauge invariance, for a comprehensive and 
detailed analysis we refer to Q: a (often neglected) term generating a consistent suppression of the decay amplitude 
is present, as will be discussed in Sections II and III . The omission of this term generates an overestimate of the 
two-photon decay rate by at least a factor of 2.25. Considering the relevance of this process related to the nature of the 
a meson, and in more general of scalar mesons (see for instance [l7l | and Refs. therein), we consider it as important 
C**** . to stress this point for future considerations about the interpretation of the enigmatic cr-resonance. 

The triangle quark loop diagram of Fig. 1 is typical for theories with quarks as effective degrees of freedom 
0, 0, B, B, El, EH EH EH]- It is evaluated both in the framework of local and nonlocal a-nn vertices. In the local 
case the Goldberger-Treiman relation on the quark level and the linear realization of chiral symmetry allow to fix the 
corresponding a-nn coupling constant. In the nonlocal case the finite size of the cr-meson interpreted as a quarkonium 
state is described by means of a covariant vertex function. The results of local and nonlocal approaches are similar 
when M a is sufficiently below threshold set by the sum of constituent quark masses. Close to threshold care has to 
be taken in the local case since a second possible problem can arise: the explicit momentum dependence of the a-nn 
coupling constant cannot be neglected further. This is also a delicate point to be treated with attention. 

The present article is organized as follows: in the next two sections we elaborate on the formalism and the results 
for a — > 77 as based on a local and nonlocal interaction Lagrangian, respectively. In section IV we summarize and 
give our conclusions. 
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FIG. 1: Quark-loop diagram contributing to Hyy decay, where H — n or a. 
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II. LOCAL CASE 



We consider the following local (L) interaction Lagrangian 

£tt( x ) = q{x)q{x) + ^q{x)i-f 5 n(x)fq(x) (1) 

where q T = (u, d) is the quark doublet of u and d quarks with the constituent mass m q — rn u — rrid (we restrict to 
the isospin limit) to be varied between 0.25 and 0.45 GeV, cr(x) and n(x) represent the scalar-isoscalar quarkonium 
and the isotriplet pion field, respectively, g a and g v are the corresponding coupling constants (which are later related 
via symmetry and low-energy considerations). We will denote the meson masses by = M n o = 134.9766 MeV and 
M a , respectively. The latter will be varied between 0.4 and 0.8 GeV. 

The decay of H = 7r°,cr into 77 is obtained by evaluating the diagram of Fig. 1. The decay width is explicitly 
given by: 



T H ^ = ^a 2 g 2 H ^M 3 H , H = n°,a, (2) 



where gHj-y = gH^ cQh^h /{2tt 2 ) is the effective if 77 coupling constant, a is the fine structure constant, N c = 3 the 
number of colors and gn refers to the coupling constant g„ or g a entering in the interaction Lagrangian of Eq. {1} . 
The charge factors Q^o = -^(| — §) = and Q a = -^(f + §) = ^= correspond to the flavor wave functions 

7T° = \J~^(uu — dd) and er = yj^(uu + dd). Finally, the loop integrals Ih [8] corresponding to Fig. 1 are functions of 
m q and Mh, which are explicitly given by 



1 2m ( M 

I„ = 4.(m„M.) = m, / d'aS | 1- » ] ^ _ = a,c 8 i„ 2 ^ ) , 
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(4) 



Note that the only difference between I w o and Jo- is the term proportional to —4a\a2 present in the integral expression 
of la, which is generally neglected in the literature (that is, the amplitudes of ir° — » 77 and a — > 77 cannot be set 
equal to each other as done, for instance, in [H, HH). The presence of the term —Aa\a2 generates a destructive 
interference with the first term, which leads to a sizable reduction of the full amplitude I a . Quantitatively the ratio 
of amplitudes is limited by I a {m q ,x) / I^o{m qi x) < 0.667 for values of < x < 2m q in the region of applicability. 
Thus, neglecting the additional term in l a implies an overestimate of the decay rate F (T _ >77 by at least a factor of 
0.667~ 2 = 2.25, as already indicated in the Introduction. Notice that we compare the decay amplitudes I a (m q ,x) 
and I„o (m g , x) but not the corresponding decay widths: as shown below these will differ consistently because of the 
dependence on the third power of the meson mass in eq. ([2]) . 

Let us now turn to the explicit calculation of decay rates. The Goldberger-Treiman (GT) relation g v = m q \f2j ' 
with = 92.4 MeV allows to determine g w . As an outcome we obtain for the decay width r, r o_ +77 = 7.73 — 8.12 
eV for constituent quark masses in the range m q = 0.45 — 0.25 GeV, in good agreement with the experimental result 
= 7.7 ± 0.5 ± 0.5 eV [20l |. Only a very weak dependence on m q is observed. 

The linear realization of chiral symmetry implies g a = g^ [H, [l2l . [l9j . Predictions for r CT ^ 77 can then be obtained 
in dependence on the effective quark mass m q and on M a . Note that we limit the parameter space by the relation 
M a < 2m q : in fact, only when this condition is met the amplitude I a remains real and no unphysical decay of the 
sigma meson into a quark-antiquark pair is included. Furthermore, the condition g a = g^ can only be employed, if M a 
is safely below the threshold 2m q . For M a ~ 2m q the momentum dependence of g a becomes non- negligible leading 
to a value for g a smaller than the one obtained in the GT limit m q y2/ ' F^; in the next section we illustrate this point 
in the context of the nonlocal approach. 



Table 1: r cr ^ 77 in the local case for 
m q = 0.25 - 0.45 GeV at M a = 0.440 GeV 



m q (GeV) 


0.25 


0.3 


0.35 


0.4 


0.45 


r^ 77 (keV) 


0.54 


0.45 


0.41 


0.39 


0.37 
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FIG. 2: ro-_ 77 in the local case as function of M a for m q = 0.35 (dark) and m q = 0.4 (gray). 

In Table 1 we report the results for T a ^^ at a fixed pole mass of M a — 440 MeV as favored by recent theoretical 
and experimental works [lj| [2(| • The results are weakly dependent on m q and clearly point to a decay width smaller 
than 1 keV, when the sigma meson is interpreted as a quarkonium state. Note for example that the omission of 
the previously discussed term in Eq. (U) implies an overestimated decay width of r CT ^ 77 = 1.18 keV for a value of 
m q = 0.3 GeV, to be compared to the correct result of 0.49 keV reported in Table 1. 

In Fig. 2 we indicate the dependence of r <T ^ 77 on M a for values of the constituent quark mass, m q = 0.35 and 0.4 
GeV, very often used in phenomenological studies. For values of M a safely below threshold (up to 0.5 GeV) results for 
r a ^ 77 lie below 1 keV and essentially do not depend on the quark mass. For values of M a approaching threshold the 
dependence on m q becomes more pronounced, where the results for rcr^ 77 eventually grow beyond 1 keV. However, 
the local approach is no longer applicable for values of M a close to threshold, as will be evident from the discussion 
of the next section. 



III. NONLOCAL CASE 

Next we the study the sigma meson described by the nonlocal (NL) interaction Lagrangian 0] 

(x) = ^=<r(x) J dhj $(y 2 ) q(x + y/2)q(x - y/2) , (5) 

where the delocalization takes account of the extended nature of the quarkonium state by the covariant vertex function 
$(y 2 ). The (Euclidean) Fourier transform of this vertex function is taken as &{k%) — exp(— k^/A 2 ), also assuring 
UV- convergence of the model. The cutoff parameter A will be varied between 1 and 2 GeV, corresponding to an 
extension of the a of about I ~ 1/A ~ 0.5 fm. Previous studies [2l[ have shown that the precise choice of 
affects only slightly the result, as long as the function falls of sufficiently fast at the energy scale set by A. The 
coupling g G is determined by the so-called compositeness condition Z a = 1 — Yl' a (M 2 ) = pLM I22I] . where T.' a is the 
derivative of the er-meson mass operator given by 

^{p 2 )^-glN c J ^-^ 2 {-k 2 )tr[S q {k+p/2)S q {k~p/2)] , (6) 

where S q (k) = (m g — ft) _1 is the quark propagator. Note, the compositeness condition is equivalent to the hadron wave 
function normalization condition in quantum field approaches based on the solution of the Bethc-Salpeter/Faddeev 
equation (23|. At this level it is clear that g a is a function of M a . In Fig. 3 we give the dependence of g a (M 2 ) at 
m q = 0.35 GeV for cut-off values of A = 1 and 2 GeV, respectively, and indicate the local limit with g G — m q V2/ F^. 
For low values of M a the coupling g a (M 2 ) is a slowly varying function, values of which for A = 1 — 2 GeV also include 
the GT limit. However, for values of M a approaching threshold g a (M 2 ) decreases below the local result (see details 
inRef. [13). 

We turn to the a — > 77 decay amplitude, where a similar suppression is found. Due to the presence of the vertex 
function $(y 2 ) inclusion of the electromagnetic interaction is achieved by gauging the nonlocal interaction Lagrangian 
([5]): in addition to the photon-quark coupling, already present in the local case, in leading order a new vertex arises, 
where the photon couples directly to the CT77 interaction vertex, see Q for details. In particular, in addition to the 
triangle diagram of Fig. 1 we have additional diagrams (see Fig. 5 in Ref. Q) to fully guarantee gauge invariance of 
the transition amplitude. In practice it is convenient to split the contribution of each diagram into a part which is 
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FIG. 3: M CT -dependence of the coupling g a {M%) at m q — 0.35 for cut-off values of A = 1 (dark) and 2 GeV (gray). The dashed 
line corresponds to the GT limit. 



gauge invariant and one which is not. The remaining terms, which are not gauge invariant, cancel each other in total 
and in the further calculation one should only proceed with the gauge invariant terms of the separate diagrams. It 
was shown [8[ , that the by far dominant contribution comes from the gauge invariant part of the triangle diagram of 
Fig. 1. The gauge invariant parts of the other diagrams are strongly suppressed (see discussion in Refs. [1, Q). 

Following [g, |9( the contribution of the gauge-invariant part of the triangle diagram to the two-photon decay width 
is given by: 



r(2) _ 



d A k 



T - /(I) + r(2) 

1 



ir 2 i (m 2 , -pf)(m| -$){m\ - p%) ' 



d A k 



4 32 

M k2 - m {kqi)[kq2) 



TT 2 i {m 2 q -p 2 ){m 2 q -pl){m 2 q -pl) 



(7) 
(8) 

(9) 



(2) 

where qi and qi are the photon momenta and p\ = k + q\, p2 = k, P3 = k — q2, q = (pi +ps)/2. The term 
contributes with opposite sign relative to 1^ leading to destructive interference. In the local limit, i.e. A — > oo, 

reduces to the term proportional to — 4a±a2 in ([!]). Note that in the pion case only a term analogous to 
contributes. 

In Fig. 4 we report the results for r CT ^ 77 in the nonlocal case as function of M a for m q = 0.35 GeV, taking values of 
A = 1 and 2 GeV. We also indicate the previous local result. While for small M a both approaches, local and nonlocal, 
agree, for increasing M a the nonlocal approach delivers smaller decay rates than the local counterpart, because of 
the threshold effects described above. The nonlocal results depend very weakly on the value A, implying that the 
numerical values for r,j^ 77 are hardly model dependent. In Table 2 we summarize our results for the two pole masses, 
M a = 0.44 and 0.6 GeV, choosing different values of m q both for A = 1 GeV and, in parenthesis, for A = 2 GeV. 




FIG. 4: T CT ^ 77 in the nonlocal case as function of M a for A — 1 GeV (dark) and 2 GeV (gray). The quark mass is set to 
m q — 0.35. The upper dashed line corresponds to the local limit evaluated in Section II. 
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Table 2: r CT ^ 77 in the nonlocal case for m q = 0.31 - 0.45 GeV, A = 1(2) GeV at M„ = 0.44, 0.6 GeV. 



m q (GeV) 


0.31 


0.35 


0.40 


0.45 


r CT ^ 77 (keV) 

at M a = 0.44 GeV 


0.238 
(0.196) 


0.192 
(0.159) 


0.152 
(0.127) 


0.124 
(0.105) 


T,^ 77 (keV) 

at M a = 0.6 GeV 


0.529 
(0.512) 


0.458 
(0.415) 


0.361 
(0.327) 


0.294 
(0.267) 



The decay widths decrease slowly for increasing quark mass while the dependence on the cutoff is very weak. The 
numerical analysis shows that 

r ff ^ 77 < 1 keV for M a < 0.7 - 0.8 GeV . (10) 

(2) 

Again, inclusion of the term ' of Eq. © is crucial to obtain these small decay widths. For instance, omission of 
this term leads to the incorrect result of r tT ^ 77 = 1.9 keV for values of m q = 0.35 GeV, A = 1 GeV and M a = 0.6 
GeV, which is almost a factor 4 larger than the correct result of 0.458 keV given in Table 2. 

In [H, using a Coulomb- like potential, the following expression relating the two-photon decay widths of tensor and 
scalar states has been derived 

r CT ^ 27 (o ++ ) = fc( ^[°^ ) rw 27 (2++) (ii) 

where m — 3. The coefficient k is 15/4 in a non-relativistic calculation, but becomes smaller (fc ~ 2) when considering 
relativistic corrections. Choosing as input M N {2++) = 1.275 GeV and r 7 j„_ >27 (2 ++ ) = 2.60 ± 0.24 keV, Eq. (Till) 
results with k — 2 in values of Tn n ^ 2l {^ ++ ) ~ 0.21 and 0.54 keV for M a = 0.44 and 0.6 GeV, respectively. These 
results are close to the corresponding numbers of Table 2. As discussed in [25j different values of the parameter m 
are obtained for different forms of the quark-antiquark potential: for instance, m = —1/3 corresponds to a linear 
potential. Then, in [25| the value m — in Eq. (| 11 1) is considered and in Ref. @ a value m = 0.3 — 1, leading to a 
larger decay width, is discussed. Here notice that our result for a light quarkonium is rather in agreement with the 
choice m = 3 and with k ~ 2, see also the model in Q where an even smaller value of k is obtained. 

Notice that we have only considered sigma masses below the constituent quark mass threshold with M a < 2m q 
and masses m q in the range of 0.25 to 0.45 GeV. In order to go beyond this limit one should (i) cither increase the 
constituent quark mass as done in Ref. Q where the 77 decays of the scalars between 1 and 1.8 GeV have been 
investigated or (ii) use more general quark propagators which include or mimic confinement. A drawback of these 
extensions is that the results reached contain a stronger model dependence, thus we do not consider these options 
here. 



IV. CONCLUSIONS 

In this work we use the formalism developed in Ref. [8j to study the decay of a scalar quarkonium state into two 
photons focusing in particular on a technical caveat of this process: a term, not present in the usual ir° — > 77 transition 
amplitude, is responsible for a sizable suppression of the a — > 77 decay rate. We considered the process a — > 77 
in the quarkonium picture both for local and nonlocal approaches. In particular the nonlocal approach allows for a 
realistic treatment of the finite size effects of the cr-meson. Similar results are obtained in both cases for masses M a 
well below the 2m q threshold. Closer to threshold the two-photon decay width in the local case should be taken with 
great care, since the momentum dependence of the coupling constant is not properly taken into account. Only the 
nonlocal result, including finite size effects, is reliable and considerably smaller than for the local case. 

Our final result r cr ^ 77 < 1 keV is smaller the results of dispersive analysis of reaction 77 — > 7r 7r° done in 
Refs. HE|- Note, that the framework developed in Ref. [26| was based on approach of Ref. The result of Ref. [f| 
evaluated at the M a = 441 MeV is r CT ^ 77 = 4.1 ± 0.3 keV, while the result of Ref. [26j is around a 40% smaller than 
that in Ref. @ , mainly due to a smaller aim coupling. 

When discussing our results it is important to stress that two aspects have not been considered. The first one is 
the possible role of pion loops. Note, that we consider a scenario where the a meson is a pure qq Fock state and, 
therefore, the a couples directly to its constituents - quarks. The coupling with other mesons (e.g. pions) goes via 
quark loops (a direct coupling of the a to pions is not present). Inclusion in a such picture of pion loops generating 
(7 — > 77 transition can occur as in Fig. 5: the corresponding amplitude is suppressed of a factor 1/N C . Our framework 
is restricted to the one-loop approximation and to the dominant term(s) in the 1/N C expansion. However, being in 
Nature N c = 3 an explicit calculation of the next-to-leading order would surely be helpful to quantify its contribution 
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FIG. 5: Pion-loop contribution to the two-photon decay within our model(s). Thik lines refer to mesons, thin lines to quarks. 

but goes beyond the scope of present paper and is left as outlook. Notice that, if we propose that the a meson is not 
pure qq state and there is also two-pion component contribution to the a meson Fock state, then we should include 
both possible intermediate states qq and 2n contributing to the two-photon transition of the a. We plan to study the 
second scenario - a being mixture of qq and 2ir - in future. 

The second aspect is related to the inclusion of finite-width effects in the evaluation of the full 77-width of the 
sigma resonance. A careful description of this point would require the precise knowledge of the propagator of the 
sigma meson dressed by pion clouds: in such a way a definition of the spectral function allowing to integrate over the 
whole mass range up to 1.27 GeV is possible [24j . Thus, also this aspect is related to pion loops and is not performed 
here. However, using trial distributions such as Breit-Wigner one and the generalized form of Ref. [24[ and varying 
the mass and the width an increase of few percent is observed. For instance, be the decay rate r CT ^ 77 = 0.458 keV 
at M„ = 0.6 GeV as in the second column of Table 2: considering a 500 MeV wide Breit-Wigner distribution the 
integrated width T a ^ 11 reads 0.66 keV. While such effects are surely important in a precision study of the two-photon 
decay width they does not change the qualitative outcome of the present paper. The corresponding integrated signals 
decay width(s) reported by [2(| are T (7 _ >77 = 3.8 ± 1.5 keV and 5.4 ± 2.3 keV, values which are not accepted as average 
or fit. Note that a confirmation of a large experimental value, contrary to usual belief, does not favor a quarkonium 
interpretation of the sigma meson. As noted in the PDG2000 [27|, the large value for r CT ^ 77 could arise from an 
additional contribution of the broad /o(1370). A clear experimental determination of the two-photon decay width 
would certainly help in clarifying the discussion related to the nature of the er-meson. 
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